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Abstract 

We calculate the branching ratio of a bottom hadron decaying into a charmed hadron and a 
charged p meson within the QCD factorization approach. We consider the effect of the finite width 
correction of the p meson. Our numerical calculation shows an obvious correction because of this 
effect. We find that the finite width effect of the p meson reduces the branching ratios by about 
9% to 11% for bottom meson decay channels: — D^p^, — )• p^, and B'^ — )• p^, and 
increases the branching ratio by about 10% for — )• A^p^ . 
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I. INTRODUCTION 



The bottom quark, as one of the third generation of quarks, was very interesting even 
before its discovery by Fermilab in 1977 [1]. The existence of a third generation of quarks is 
necessary for the complex phase to exist in the Cabibbo-Kobayashi-Maskawa (CKM) matrix, 
which is the source of CP violation within Standard Model j2|. Because of this, studies of 
the physical processes involving the bottom quark are very important for testing Standard 
Model and for finding New Physics. Non-leptonic weak decays of bottom hadrons supply 
a rich platform to dig the properties of the bottom quark. The technique of Operator 
Production Expansion allows us to calculate decay amplitudes of heavy hadrons via an 
effective Hamiltonian jsj. In the calculations, one always confront hadronic matrix elements, 
which contain non-perturbative effects and are hard to handle. However, with the aid of 
Heavy Quark Effective Theory [j] , one can simplify the hadronic matrix elements via several 
approaches, such as perturbative QCD js], QCD factorization js-Q], 



and Soft Collinear 



Effective Theory [lo|4l3|. 



The charmed decays of bottom mesons have been studied extensively in experiments 
for quite a long time. For example, CLEO reported that the decay branching ratio for 
B° D+p- and B+ D^p+ are 0.0135 ± 0.0012 ± 0.0015 and 0.0077 ± 0.0013 ± 0.0002, 
respectively [141] . 

The properties of bottom baryons are less clear to us than bottom mesons mainly due 
to its rare production rates compared with bottom mesons. However, the decays of A;, have 
been becoming more and more clear. The first channel of nonleptonic decay of yet has 
been measured is A^ — )■ Anir 15 1 . Besides, CDF CoUaberation also measured two charmless 

n 

decay channels of A^: A;, — > pir and Af, — )■ pK [16|]. The running of Large Hadron Collider 
(LHC) gives us a very good opportunity to improve the experimental sensitivity. In the 
near future, LHCb will collect enough events to observe more decay channels of A;,. So it is 
urgent for us to study decays of heavy baryons such as A^ theoretically. 

In this paper, we will focus on one specific kind of decays of the bottom hadrons, in which 
a bottom hadron decays into a charmed hadron and a charged p meson. Since the p meson 
decays into two pions strongly, we will take into account its finite width effect. To our latest 
knowledge, in most of theoretical considerations for such kind of decays usually a zero width 
limit for the p meson is used. This leads to less accurate results because the decay width of 



the p meson is not small. 

This paper is organized as follows. In Sec. HIl we briefly give the formalism for the 
decay processes we deal with. In Sec. IIIIl we show the numerical results and make some 
discussions. The conclusion is presented in Sec. HVl 



II. FORMALISM 



Our start point is the effective Hamiltonian for the weak decay h — ?■ udc 

^eff = ^^udVcb [CMQM + C2(/i)Q2(/i)] + H.C., 



(1) 



where Gi? is the Fermi constant, ci and C2 are Wilson coefficients at the scale /i (which is 
(9(mb)), Vud and Vcb are CKM matrix elements, Qi and Q2 are four quark operators which 
take the form 

where u, d, c, and b represent quark field operators, the superscripts i and j are color indices. 
The Wilson coefficients ci(/i) and C2{p) take the form 
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In the above equations mw is the mass of W boson, k is a scheme dependent parame- 
ter, which equals 0(=f4) for Naive Dimension Regularization ('t Hooft-Veltman) scheme jsl. 



Noticing z±{fi) is independent of renomalization scheme, we can define scheme independent 
Wilson coefficients as 



(6) 
(7) 



We will first consider the bottom hadron decay A° — )■ A^p" . This process was considered 
in the naive factorization approach long time ago Here, we will use the QCD factoriza- 
tion approach which was originally constructed for B meson decays w^- When calculating 
the decay width for A° — )• A+p^, after summing over polarizations of A;, and Ac, one has to 
typically deal with 



m^tr 



(8) 



where e is the polarization vector of the p meson, rup {rn\^^^^) is the mass of p {Ab(c)), Vk^^^) 
is the momentum of A{,(c)- For longitudinally polarized p meson (A = 0), one can replace 



(A=0) 



by iq {q is the momentum of the p meson), leading to 



m^tr 



While for transversely polarized p meson (A = ±), one has 



m^tr 



(A=±), 



2 2 



(9) 



(10) 



Thus the production of a longitudinally polarized p meson dominants. When we consider 
QCD corrections to the decay A° — > A+p~, only the vertex corrections contribute in the 
heavy quark limit |l^. The decay amplitude for a longitudinally polarized p~ in the final 
state is then 

Gp 



A 



where 



aij = ci(mfe) + 



V2 



KdVcbfpiAclc^iaiv - aiAl5)b\Ab), 



dx^p{x)Fj{x, z) 



(11) 



(12) 



here j = V, A, z = mc/mi,, $p(x) is the light cone distribution amplitude of the p meson. 



The form of Fj (z) can be found in Ref . 



18|. 



In the heavy quark limit, the hadronic matrix elements (A^ 1 07^,(7, 6 l A^), which corre- 



sponds to the weak transition from A5 to A^, can be parameterized as 

(Ac|c7^(75)6|Ab) = C {UJ{S)) UASpAj-ff,{-f5)uA,{pA,) 
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(13) 



where ( is the Isgur-Wise function 
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ml + ml 



22|, with 

s 



2m m ' 

where s = {q = pa^ — PaJ- Now the nonperturbative effects of strong interaction are 
fully described by the decay constant of the p meson and the Isgur-Wise function. Then the 
decay width for A° ^tP~ 



^^2 

32^ 



VuM'f.ml (1 - z'Y IC {coimD) \' {\a^y\' + \a^^\') . 



(15) 



As will be seen later, the above expression for the decay width of A° — t- A^p^ actually 
corresponds to the situation where the decay width of the p~ meson goes to zero. However, 
the vector meson p~ has a relatively large decay width (comparing with its mass) because 
it decays rapidly into two pions. In addition, the mass of the pion is not small compared 
with the mass of p~ meson. As a result, we should take the decay width of the p meson into 
account. This means we have to deal with the decay chain A^ — )■ A+p^ A+tt^tt^. 

The decay amplitude of the decay chain A° — )■ A+p^ A+tt^tt^ can be parameterized as 



M 



A0^A+p-(^7r07r-) 



where Sp = s — m'l{s) + iy/sTp{s), and 



AO^A+p- 



-^P-^ttOtt-' 
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^KdVcbVsfp{Ac\c7''iaw - aiA^5)b\Ab), 



(17) 



with s being the invariant mass square of the pion pair, ^^(s) the running mass squared 
of the p meson, T p{s) the decay width of the p meson, gpT,.„ the effective coupling constant 
among p and two pions. The s-dependence of the running mass squared m'i{s) take the form 
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for s > 4m^. The decay width of the p meson depends on s via [23|, |24| 

/ o \ 3/2 



where Tp is the nominal decay width of the p meson, which can be read directly from 
the particle property list in Particle Data Group. In general, taking into account the s- 
dependence of the effective coupling constant gpT^, ^p{s) should take the from 

3/2 
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where gp.„T^ = gpT,T,{rn?p). Through out this paper, we simply assume that (7^^,^ is independent 
of s. The coupling constant gpTr-w, corresponding to p decaying into two pions, takes the form 

AMw?pTp^2-. 487r 



2 

9 pTTTT 



{m?p — Ami) 



3/2 
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p-)-27r 



(22) 



Because the p meson decays dominantly into two pions, we will simply assume that rp_^27r 

The differential decay width for — )■ A+p^ — )■ AcTt^tt^ is then 

2 



~ (2^0^32^^^ 5Z |^-^AO-.A+p-(-..o.-) 



dsds', 



(23) 



spins 

with s' being the invariant mass square of Ac and vr^. Integrating over s', we get (in the 



heavy quark limit) 

dr 

dJ " 



ds'- 
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^/sTp{s)/7i 



dsds' [s — m2(s)]2 + sT'j^{s) 
where the lower and upper bounds of the integral are 
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For the total decay width, we have to integrate over s. The bounds for this integral are 4m^ 



and (m\ 



respectively. Although the use of QCD factorization is not appropriate 



when s is close to the upper bound of the integral, it can be seen from Eq. fl24|) that the 
differential decay width is very small when s is in this region, so the contribution of this 
part is neghgible. As a result, the decay width for A° — )■ A^p~ — )■ A^vr+vr" can be expressed 

as 



(26) 



where 



'6 



C (c^K)) 



(27) 



Note that if we let the decay width of the p meson goes to zero, we will have the Breit-Wigner 
form as 

lim "^^^'i'^^""^,, , = 6{s - ml). (28) 

So we find the decay width becomes just r° in this limit. 

In deriving Eq. 0271) . we have applied the heavy quark limit except the bound of integral 
over s. However, in order to improve numerical accuracy, we will not use this equation in 
the following. Instead, we will use 



d.d.-^, . ' .^J^X-^ Jh ^ (29) 



167r2 (m|^ - m|J3 \^{m^^) J [s - m2(s)]2 + sTl{s) 

■ [{ml^ + mljis' + s") - As's" - 2m^(mi^ + m^J + Amt] . 

Similarly, when dealing with bottom meson decays: i? —t- -Dp, if we take the decay width 
of the p meson to be zero, we find that the decay width is 



2 ^_^3 

V:,V^MDp)\'m^{K))?- (30) 



^^^P ~327r Am%mD 



However, if we take into account the decay width of the p meson, again have to deal with 
the decay chain B — )• Dp — > Dim. The differential decay width is then 



-(27r)3 32m|'^^^^''(^^^)' 



1 (^B + ^d)^ ^2 IT/ T/ |2 £2\ ( N\l2 / li t\23p-KTT 



(31) 



Similarly, the decay width is modified to 



B-)-Z)p) 

where 



^B^Dp — Rb^Dp X ^B^Dp} 

(32) 



Rb^Dp I dsds'^g^2 _ ^2^)3 (/(U)) [s - m2(s)P + sTlis' ' ^^^^ 



^ Note that in this equation, we neglect a term proportional to (|aiyp — |aiAp)/(|aiy P + |aiAp) because 
it is numerically very small 10"^). 
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III. NUMERICAL RESULTS 



Our input parameters are from Paticle Data Group 2012 25|. With the aid of the 
renormahzation group equation for the running couphng constant, we can get Aq^j-,(MS) = 
231 ± 9 MeV at 2-loop order, and Agij3(MS) = 213 ± 8 MeV at 3-loop order, where the 
uncertainties come from as{mz) and mz- Note that at 4-loop order Aq^j-,(MS) = 213 ± 
8 MeV {25!. It can be seen that the difference between 3-loop order and 4-loop order results 
is small. In fact, we can use either the 2-loop order result or the 3- or 4-loop order result to 
get the numerical values of the Wilson coefficients, the difference among the obtained results 
are negligible. Our numerical results for the scheme independent Wilson coefficients are 

ci(mb) = 1.146 ± 0.003, C2{mh) = -0.312 t 0.005, (34) 

where the uncertainties come mainly from the uncertainties of the QCD scale and rrib. The 
numerical results for the coefficients aw and aiA are then 

Re(aiy) = 1.057(2) - 0.0064a^(/i) + 0.0029a^(/i), 
Im(aiy) = 0.0223(6) + 0.032a?(/i) - 0.0015a^(/i), 

(35) 

Re(aiA) = 1.056(2) - 0.0091a?(/i) + 0.0005a^(/i), 
Im(aiA) = 0.0146(4) + 0.028<(/x) - 0.0021a^(/i). 

where and ctg are Gegenbauer moments for the p~ meson jsl. We only keep the uncertainty 
of the leading term in the Gegenbauer expansion of the p~ meson wave function, which comes 
from the uncertainties of the QCD scale and the masses of bottom and charm quarks. We 
use the following Gegenbauer moments for the p~ meson 8|: 



0.3 ±0.3, a^ = 0.1±0.3. (36) 



We use the decay channel p° — )■ e+e to extract the decay constant of the p meson, fp. 
The decay width for this channel can be expressed as 

r,o^...- = (37) 

which leads to fp = 219 ± 1 MeV. 

For the Isgur-Wise functions, it can be parameterized as 

^co) = (38) 



where is the slope parameter, and ^(1) = 1 in the heavy quark hmit. One can fit these 
parameters with certain experimental data. For B ^ D transition, according to Heavy 
Flavour Average Group (HFAG), these parameters are fitted as ^(1) x |V"cb| = (42.64 ± 



1.53) X 10"^ p2 = 1.186 ± 0.054, with a correlation of 0.829 j26|. For transition, 
the uncertainties of the fitted parameters are quite large. For exarnple, the fitted value for 
the slope parameter is = 1.59 ± 1.10 by DELPHI Collaboration 27|. 

We can also adopt some models for the Isgur-Wise function. The first model we will use 
is the MIT bag model which was proposed by Sadzikowiski and Zalewski 28| . In this model, 
the Isgur-Wise function takes the form 



5M=^^J , (39) 

where parameters (a, b) take the values (2, 0.6), (2.7, 0.6) and (3.5,1.2) for B ^ D, Bs Dg 
and Afc — > Ac transitions, respectively. For Aj, — > Ac transition, we use a so-called soliton 
model proposed by Jenkins, Manohar, and Wise 29|, in which C,{ui) = 0.99e"^'^'''^^^\ 

The numerical results are shown is Table [H We can see that the finite width effect of the 
p meson reduces all the branching ratios of the three decay channels of bottom mesons by 
about 10%. On the other hand, the finite width effect of the p meson enlarges the branching 
ratio for A^ — t- A+p^ by about 10%. The ratio R is insensitive to Isgur-Wise function 
models. The branching ratios depend on Isgur-Wise models more sensitively. We can also 
see from Table |T]that for B^ — t- D^p^ and Bg — ?► Djp^ the finite width effect of the p meson 
makes the branching ratios closer to the center value of the experimental data while for 
B^ — i- D^p^ the situation becomes worse. In our calculations we worked in the heavy quark 
limit while dealing with the warm transition matrix elements. The order-l/mfe corrections 
may numerically lead to about 10% change for the decay widths. This is beyond the scope 
of our present work. 



IV. CONCLUSION 



In this paper, we calculate the decay widths of the decay channels Af, — )■ AcP and B ^ Dp. 
We consider the QCD corrections to these decay channels. In the heavy quark limit, only 
vertex corrections contribute, and the production of the longitudinally polarized p meson 
dominants. In the heavy quark limit, the noperturbative effects of strong interactions are 



9 



fully described by the decay constant of the p meson and a form factor corresponding to the 
weak transition A;, — t- Ac or i? — t- 

Since the p meson decays dominantly into two pions through strong interaction, it has 
a very broad decay width comparing with the mass of the p meson. Also, the mass of the 
pion is comparable with the mass of the p meson. Even in the heavy quark limit, the two 
ratios Tp/nip and mT^/rrip do not go to zero. As a result, we do not neglect the mass of the 
pion meson in the calculation. 

The correction of the finite decay width effect is described by the parameter R in this 
paper, which is the ratio of the decay width with and without the finite width effect of the 
p meson. We find that R = 0.89 — 0.91 for B — ^ Dp, indicating a decrease of the branching 
ratio by about 9% to 11%. For A;, — )■ A^p channel, we find that the branching ratio is 
increased by about 10%. 
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TABLE L Branching ratios for different decay modes. IW stands for Isgur-Wise function. 



Modes 


IW R 


B = RB°{%) 


exp. B (%) 




MIT bag 1.11 1.02 


1.13 




Ab Acp 


soliton 1.07 1.74 


1.86 






DELPHI 1.10(7) 0.57-3.44 


0.67-3.55 




B+ 'Wp+ 


MIT bag 0.89 1.13 


1.00 


T.34±0.18 [l£ 


HFAG 0.90(1) 1.05(6) 


0.95(5) 


D-p+ 


MIT bag 0.89 1.04 


0.92 


■0.78 ±0.13 [14] 


HFAG 0.90(1) 0.97(5) 


0.88(5) 




MIT bag 0.89 0.78 


0.70 


0.74 ±0.14 [30] 
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